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Symmetry-protected topological phases (SPT) are short-range entangled gapped states protected 
by global symmetry. Nontrivial SPT phases cannot be adiabatically connected to the trivial disor¬ 
dered state (or atomic insulator) as long as certain global symmetry G is unbroken. At low energies, 
most of two-dimensional SPTs with Abelian symmetry can be described by topological quantum 
field theory (TQFT) of multi-component Chern-Simons type. However, in contrast to the fractional 
quantum Hall effect where TQFT can give rise to interesting bulk anyons, TQFT for SPTs only sup¬ 
ports trivial bulk excitations. The essential question in TQFT descriptions for SPTs is to understand 
how the global symmetry is implemented in the partition function. In this paper, we systematically 
study TQFT of three-dimensional SPTs with unitary Abelian symmetry (e.g., x Zjv 2 X • • •). In 
addition to the usual multi-component BF topological term at level-1, we find that there are new 
topological terms with quantized coefficients (e.g., a 1 A a 2 A da 2 and a 1 A a 2 A a 3 A a 4 ) in TQFT 
actions, where a 1 , a 2 ,-- are 1-form U(l) gauge fields. These additional topological terms cannot 
be adiabatically turned off as long as G is unbroken. By investigating symmetry transformations 
for the TQFT partition function, we end up with the classification of SPTs that is consistent with 
the well-known group cohomology approach. We also discuss how to gauge the global symmetry 
and possible TQFT descriptions of Dijkgraaf-Witten gauge theory. 


I. INTRODUCTION 
A. Background and motivation 

At low energies, universal properties of quantum many- 
body systems may be governed by continuum quantum 
field theory. For symmetry-breaking phases, Ginzburg- 
Landau field theory with symmetry-breaking potentials 
is utilized, where the concept of local order parameters is 
introduced. For topological phases of quantum matter, 
such as fractional quantum Hall effects (FQHE), local 
order parameters vanish and thereby fail to character¬ 
ize the ground state. However, such phases can be de¬ 
scribed by topological quantum field theory (TQFT) 1,2 . 
As a renormalization fixed point theory, TQFT partition 
functions are usually formulated in the Euclidean path- 
integral formalism where the classical action is purely 
imaginary and all correlation functions of local opera¬ 
tors are invariant under diffeomorphism of the spacetime 
manifold. In this sense, TQFTs are expected to effi¬ 
ciently capture the global phenomena that are insensi¬ 
tive to local energetic details. For example, the (2+1)- 
dimensional Chern-Simons field theory was introduced to 
describe FQHE 1 12 . It reveals the key mechanism of flux- 
charge attachment for FQHE and successfully predicts 
the anyon statistics and quantized Hall conductance 2 . 
Very recently, there is intensive ongoing research focus¬ 
ing on the interplay between symmetry and topology in 
strongly correlated systems. For example, the so-called 
“symmetry-protected topological phases” (SPT) in inter¬ 
acting bosonic / spin systems 13 35,37 64 have been draw¬ 
ing much attention. In contrast to FQHE where the 
bulk is fractionalized and supports topologically anyonic 


quasiparticles, the bulk of bosonic SPT phase is non- 
fractionalized and only supports trivial bosonic quasi¬ 
particles. Nevertheless, distinct SPT phases can still be 
identified as long as a symmetry group, say, G is unbro¬ 
ken. In other words, two distinct SPT phases with G 
cannot be adiabatically connected to each other unless 
G symmetry is broken 14 . 

Since TQFT approach has been successfully applied 
in topological phases of quantum matter (e.g., FQHE), 
a natural question would be: Can we also use TQFT 
to describe universal properties of SPT phases that, by 
definition, do not support fractionalized bulk excitations? 
If so, we expect that such TQFTs must have two key 
properties: 

1. Bulk excitations are non-fractionalized. In other 
words, they have a unique ground state on any 
closed manifold and there do not exist nontrivial 
topological sectors. 

2. TQFTs of distinct SPT phases are topologically 
distinguishable from each other if and only if cer¬ 
tain unbroken symmetry G is considered. It indi¬ 
cates that nontrivial properties of TQFT are topo¬ 
logically robust only when G symmetry is unbro¬ 
ken. 

In Ref. 18, it was proposed to use multicomponent 
Abelian Chern-Simons theory to describe 2D SPT phases 
with Abelian symmetry. The Chern-Simons theory ap¬ 
plied in SPT phases must have a unique ground state 
on a torus. SPT phases are classified by distinct ways 
of anomalous symmetry transformations implemented on 
the ID boundary. In Ref. 25, a complete TQFT descrip¬ 
tion for 2D SPT phases with Abelian symmetry has been 
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FIG. 1. (Color online) Bosonic wormhole effect as a physical 
mechanism of SPT phases with J^Za^ symmetry, (a) A 
ZjVi symmetry domain wall “D(ZjVi)” and aZjv 2 symmetry 
domain wall U V(Z n 2 )” intersect along a closed loop “C” (only 
a segment is shown due to the limited space), (b) shows the 
bosonic wormhole effect induced by intersections of symmetry 
domain walls. Each intercepting line forms a closed loop in the 
bulk, but some of them may form open strings and end at two 
2D boundaries between which the bulk is sandwiched. Along 
such an open string C (the line with arrow), Zjv 2 symmetry 
charge Q is pumped from the endpoint on the lower boundary 
to the endpoint on the upper boundary. See main text for 
more detailed explanation. 


TABLE I. A brief summary of irreducible 3D SPT phases with 
unitary Abelian symmetry. By “irreducible”, we means that 
all subgroups of symmetry group play nontrivial roles in pro¬ 
tecting the nontrivial SPT phases, a 1 and b 1 are 1-form and 
2-form U(l) gauge fields, respectively. denote the 

corresponding classifications, where Nij... are greatest com¬ 
mon divisors of Ni,Nj , • • •. SPT phases with either Zn or 
U(l) fc or Z n x U(l) fc are trivial and not included below. All 
other SPT’s with unitary Abelian group symmetries can be 
obtained directly by using this table. For example, with arbi¬ 
trary unitary finite Abelian group G — Z n x x Z n 2 * * *, it is suf¬ 
ficient to use this table to reproduce all SPT phases classified 
by the group cohomology, since 'H 4 (Zn 1 x Zn 2 * * * , U(l)) = 
Ll/< j(ZjVjj ) 2 x IIj< j<k(Zn ijk ) 2 x U 1 <J<K<LZNUKL • 


Symmetry Topological Quantum Field Theory and 
Group Classification 

Z jVi x Zjv 2 

2F / Yi b 1 A da 1 + ipi f a 1 A a 2 A da 2 (Zjv 12 ); 

^ I Y^t b 1 A da 1 + ip2 f a 2 A a 1 A da 1 (Zn 12 ) 

Z N 1 x 
Zjv 2 x Zjv 3 

2F I Y"i b 1 A da 1 + ipi f a 1 A a 2 A da 3 (Zjv 123 ); 

^ I Y 3 r b 1 A da 1 + ip2 f a 2 A a 3 A da 1 (Ziv 123 ) 

ri/ zjvj 

2 ^ / Yi b 1 A da 1 +ip f a 1 A a 2 A a 3 A a 4 (Z n 1234 ) 

Zn 1 x 

Z n 2 x U(l) 

2F / Yi b 1 A da 1 + ip f a 1 A a 2 A da 3 (Zjv 12 ) 


achieved. Nevertheless, a complete TQFT description of 
3D SPT phases is more challenging and still not con¬ 
clusive even for Abelian symmetry groups. In Ref. 55, 
bosonic topological insulators (BTI) as specific cases of 
3D SPT phase are studied and the corresponding TQFTs 
are derived via the vortex-line condensation scenario, 
where the cosmological constant-type term “6A6” is used. 
Ref. 53 proposes to use topological BF field theory to de¬ 
scribe some SPT phases in the presence of U(l) charge 
conservation symmetry. But except these special cases, 
we actually know very little about the TQFT description 
of 3D SPT, especially with unitary discrete symmetry. 
Indeed, there have already been some interesting propos¬ 
als of 3D SPT materials, such as “topological paramag¬ 
netism” in 3D integer-spin Mott insulators 65 . We hope 
that a full establishment of TQFT in every spatial di¬ 
mension will help us to systematically identify low-energy 
properties of bosonic SPT phases and realize them in a 
lab in future. 


alizing GL actions, we obtain TQFT partition functions 
(Sec. III). We quantitatively study how symmetry trans¬ 
formations are realized in TQFTs, which gives rise to the 
same SPT classification results that were previously ob¬ 
tained from group cohomology approach 14 . While SPT 
phases with JI/ ^iVj are studied in details, other SPT 
phases can be understood in a similar manner (Sec. IV), 
which are summarized in Table I. As such, TQFTs in 
this paper may be regarded as continuum field theory of 
group cohomological lattice models. We also connect our 
TQFT results to the “decorated domain wall picture” 40 
as well as the interesting ‘3-loop statistics’ phenomena 62 
(Sec. IV). A brief summary of the paper including some 
future directions is given in Sec. V. Several technical de¬ 
tails are collected in Appendices. 

II. SPT AS A QUANTUM DISORDERED 
PHASE: GINZBURG-LANDAU ACTION 
APPROACH 


B. Main results and outline 

In this paper, we attempt to establish a TQFT de¬ 
scription of 3D SPT phases with unitary Abelian symme¬ 
try (see Table I) . Instead of the abstract group cohomol¬ 
ogy approach 14,17 , the strategy we adopt in this paper is 
based on Ginzburg-Landau (GL) type actions with global 
symmetry (Sec. II). Such an approach is more physical 
and accessible. The physical mechanism for realizing 
SPT phases is pictorially illustrated in Fig. 1, where the 
“bosonic wormhole effect” is introduced. Then, by du- 


A. A review of quantum disordered phase 

We start with a brief review of the quantum field 
theory description for quantum disordered phase. Let 
us consider the following partition function of four¬ 
dimensional classical XY model in Villain form: 



where G Z is an integer link variable defined on the 
NNN link (ij) of 4D cubic lattice, x is inverse tempera¬ 
ture. Oi G (0, 27t] is a U(l) phase angle defined on lattice 
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site i. This classical model corresponds to a quantum 
XY model in 3 + ID spacetime. There is a phase transi¬ 
tion which separates two quantum phases at x = Xc- At 
small x region near Xc where the Villain approximation 
works, we have a “quantum disordered” phase where the 
ground state respects U(l) global symmetry and the low¬ 
est excitations are gapped. Such a quantum disordered 
phase serves as our starting point for understanding SPT 
phases. 

To proceed further with field theory language, we 
rewrite the above lattice theory in terms of continuous 
spacetime variables with the following Lagrangian: 

£o = §(^0) 2 = |(a^ 8 +V«) 2 = | (duOs-arf, (2) 

where we decompose 0 as 6 = 0 S + 0 V . The 0 S variable in 
Eq. (2) is a smooth function. The singular part (or vortex 
part) 0 V of original lattice model can be incorporated 
in the newly-introduced vector variable a M . It can be 
viewed asaU(l) gauge field since the action is manifestly 
invariant under 6 —x 6 + y, —)> + d^x where x is a 
scalar gauge parameter. Apparently, such a gauge degree 
of freedom origins from the ambiguity(up to the scalar 
field x) of the decomposition 9 = (9 S + x) + (Pv ~ x). 

We note that the action Eq. (2) actually only describes 
one of local minima of the Villain form Eq. (1). Since 
9 is by definition smooth, the large gauge transforma¬ 
tion of a with singular x and f sl d\ = 27rn (with n an 
arbitrary integer) will transform the action Eq. (2) to 
another local minimum. Therefore, in order to main¬ 
tain large gauge invariance, it is required to sum over 
all local minima, which essentially reproduce the Villain 
form Eq. (1). Thus, in the path integral formulation 
Z = f D0Dae~ f d xCo , all the large gauge transforma¬ 
tions of a will be naturally included. 


B. Ginzburg-Landau action description for a 3 + ID 
bosonic SPT phase with G — ~Ln 1 x Zjv 2 symmetry 


Since there are no interesting SPT phases with Z^v or 
U(l) fc symmetry in 3 + ID, here we consider the sim¬ 
plest 3 + ID bosonic SPT phase with Abelian symmetry 
G = x Zat 2 . We begin with the following Ginzburg- 
Landau action with U(l) x U(l) symmetry. 


S Q = 



1 

4 g 2 


(A ) 2 


+ •••• (3) 


At the quadratic level of the low energy effective the¬ 
ory, we may also add a “Maxwell” term A_(/‘^) 2 where 
fixv — d^a u — and g is gauge coupling constant. 

In addition to So, we may add the following topological 
term: 

Stop = ~ip jd 4 xe ,J ' uXp (d^9l - a* )(d„0 2 - al)d\a 2 p . (4) 


Mathematically, Sx op is a wedge product of differential 
forms that is invariant under diffeomorphism. Sx op is 
the only term in (3+l)D spacetime that is topological 
and gauge-invariant with two independent 1-form gauge 
fields. Although Sx op formally has a U(l) x U(l) symme¬ 
try, below we will show that such a symmetry is anoma¬ 
lous and can not be realized in a strictly 3D system. 
However, when p takes certain quantized value, Eq.(4) 
has an anomaly-free Z^q x Z n 2 symmetry. In Sec. Ill, 
we will rigorously show that p is not only quantized but 
also takes values in a compact space. Distinct SPTs may 
be labeled by distinct p’s, which may unveil the classi¬ 
fications of SPT phases. As a matter of fact, a usual 
dimension reduction scheme may already provide us a 
simple picture about the quantization condition on p. 

For this purpose, let us consider a special 4D manifold 
A4 4 = T 2 x T 2 , where f T2 da 2 = 2tt on one of the torus 
T 2 . After compactifying over this T 2 , we end up with a 
1 + ID SPT phases described by the following topological 
term: 

Stop = ~ 27ri P Jd^xe^id^l - al)(d v 0 2 s - a 2 ). (5) 

In Ref. 25, it has been shown that if 47r 2 p = 0 mod 1, 
the above topological term has an anomalous U(l) x U(l) 
symmetry and can not be realized as a pure 1 + ID sys¬ 
tem. However, if 47r 2 p = 0 mod N ^ 2 2 with Ni 2 the 
greatest common divisor of N\ and 7V 2 , the above topo¬ 
logical term has an anomaly-free x I N2 symmetry. 
Therefore, if p = where k = 0,1, • • • Vi 2 , the orig¬ 

inal topological term Eq. (4) will describe Vi 2 different 
3 +ID SPT phases with Z^ x Z n 2 symmetry. In Sec. Ill, 
the classification will be systematically derived from a 
TQFT which is a dual theory of the above GL action. 

The physical meaning of St op can be understood via 
bosonic wormhole effect induced by intersections of sym¬ 
metry domain walls. On one hand, the wormhole effect 
of topological insulators was discussed in Ref. 69, where 
the bulk is a fermionic system. On the other hand, by 
definition, a Z^ symmetry domain wall appears as a 
2D closed manifold that separates two 3D regions with 
two distinct 6 1 = j^kj with ki E Za^. We consider 

= a 2 = 0 configurations such that Z jsf T symmetry do¬ 
main walls (/ = 1, 2) do not proliferate. As a result, there 
are many closed loops that are intersections between Z ^ 
and Zjv 2 symmetry domain walls as shown in Fig.l-(a). 
Such closed loops do not contribute a Ztv 2 symmetry 
charge Q from Stop since Q = di(dj0 1 dk0 2 )e ljk = 0. 
However, if open boundaries are considered, intersections 
may form open strings that end at two 2D boundaries 
as shown in Fig.l-(b). At these endpoints, nonzero Q 
is realized as a pumped charge due to wormhole effect. 
Each symmetry domain wall is labeled by 27 ^~ J , where 

Ki e T Ni (7 = 1,2). As such, Q=p 2 ^' = k -^ 

(k E Zat 12 ) where the quantization condition Eq. (14) is 
applied. N 12 is the greatest common divisor of N± and 
V 2 . Once a* and a 2 are turned on, the energetic cost of 
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symmetry domain walls is compensated by gauge config¬ 
urations. Thus, the above novel domain wall configura¬ 
tions proliferate leading to a symmetric ground state. 


III. TOPOLOGICAL QUANTUM FIELD 
THEORY OF SYMMETRY-PROTECTED 
TOPOLOGICAL PHASES 


A. Duality between GL theory and TQFT 


In the following, we may apply duality transformations 
to obtain a TQFT description of GL action. We start 
with the continuum action So + -Stop and derive TQFT. 
Next step is to apply Hubbard-Stratonovich transforma¬ 
tion: 

c =i J2KK - 5 X) + E %r- 

i i y 

- ipe^ xp {al - d^e\)(al - d v 0 2 s )d x a 2 p + T ^(n') 2 , 


where the vector fields n 1 are auxiliary fields in order 
to linearize the quadratic terms. After total deriva¬ 
tive terms are dropped off, integrating out 6\ leads to: 

— pa 2 d\a 2 e puXp ) = 0 which can be resolved by 
introducing a 2-form gauge field b x pv \ (b 1 ^ = — bj, ) 

nl = pald x ay" xp + . 

After replacing n p by the above identity, we end up with: 


£ = in K a l ~ d Al) + i jz a l d ^ b i P e ' luXp 


An 

■ ipe puXp a]\d v 0 2 s d\a 


2 

E 

i 


I \2 


uy 

4 g 2 


E^«) 2 ’ 


where the last term will be replaced at the final step. 

Now, further integrating out 0 2 leads to: c^(n 2 + 
pald\a 2 e puXp ) = 0 which can be resolved by introduc- 
ing a 2-form gauge field b(b 2 ^ = -6^) 

nl = -paldxay^ + ^d u b\y vXp . 

After replacing n 2 by the above identity, we end up with 
the partition function with the TQFT action: 


Stqft = /*A da Ij r ip Ja 1 A a 2 A da 2 + Sm, (6) 


where b pi/ are three 2-form U(l) gauge fields. As a side 
note, the differential form notation ‘fr 7 ’ is replaced by 
once spacetime indices are written explicitly. 


To the best of our knowledge, the term f a 1 A a 2 A da 2 
as a classical action was introduced in Ref. 73 where two 
additional constraints da 1 A da 2 = 0 and da 2 A da 2 = 0 
are imposed to recover gauge invariance 73 . And as 
an ultraviolet regulator, is given by: 

S M =J d 4 x A ^^ay " xp + ±-d v b\y ' x '>) 2 

+1 d 4 xA(_ pa i dxa y^P + l^ dvh i pe ^ P f 

Interestingly, from Sm, one can explicitly read the gauge 
transformations (9) that play a very important role in 
the remaining discussions. 


B. Exotic gauge transformations 

The first term in Eq. (6) is two copies of topological BF 
term 70 72 at level-1. As usual, all gauge fields are subject 
to the following usual Dirac quantization conditions: 

T f da'ezAf db i G z, (8) 
Jm 2 27T Jm 3 

where M 2 is any closed two-dimensional surface. As 
an Abelian gauge theory, the gauge transformations are 
given by: 

a 1 a 1 + dx 1 , b 1 ^ b 1 + dV 1 - 2npe IJ3 X J da 2 , (9) 

where X T G R and V* G 1R are independent scalars and 
vectors, respectively. Here, both 1-form d X T and 2-form 
dV 1 are closed but allowed to be non-exact if the follow¬ 
ing integers are nonzero on non-contractable manifolds. 

f dx 1 = 2nn I , f dV 1 = 2irk 1 , (10) 

Jm 1 Jm 2 

where n 1 and k 1 are six independent integers. For 
nonzero n 1 and k 1 , the corresponding gauge transforma¬ 
tions are said to be “large gauge transformations” labeled 
by the winding numbers n 1 and k 1 . 

If p = 0, b 1 are transformed in the usual definitions of 
2-form gauge fields. Remarkably, the presence of p term 
in Eq. (6) induces a p-dependent term in the gauge trans¬ 
formations. Due to such unusual gauge transformations, 
we may construct the following two gauge-invariant op¬ 
erators: 

Wilson loops: e 1 ^ 1 a 

Modified Wilson surfaces: b ' ~ i27rp eIJ3aJAda2 

where V 3 is an open volume enclosed by Ad 2 , i.e., dV 3 = 
M 2 . It means that the Wilson surface operators of b 1 and 
b 2 are not gauge invariant alone, in contrast to the usual 
definitions. In Fig. 2, these gauge invariant operators are 
shown pictorially. 
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FIG. 2. (Color online) Illustration of gauge-invariant oper¬ 
ators introduced in Sec. IIIB. (a) The Wilson loop opera¬ 
tor e z ^-M ia . (b) The operator formed by exp{z J m2 b 1 — 

i2vrp f v3 e IJ3 a J A da 2 } where M 2 — dV 3 . In (b), the cube rep¬ 
resents V 3 and its surface represents M 2 . The star symbols 
in (b) represent nonzero contributions of the “Chern-Simons 
density” f v3 e IJ3 a J A da 2 in V 3 . 


C. Global symmetry G = Tn 1 x Zjv 2 in TQFT 


In TQFT action (6), the physical meanings of the 
two form gauge field b can be understood by identify¬ 
ing its curvature as four-current of point-particles: = 

_L ^vXpq^i i n terms of differential forms, we have: 

* jI = T^db 1 , ( 11 ) 

where ★ denotes Hodge dual operation. If the symmetry 
is U(l), then, we may add a U(l) external background 
gauge potential A 1 to impose the U(l) symmetry: 


^TQFT 


S - Stqft + E 2^ / 


— / A 1 A db 1 


( 12 ) 


where constant terms are neglected unless otherwise 
specified. In our case of G = the particle cur¬ 

rent is conserved mod Nj. To impose this condition, the 
closed loop integral of the 1-form A 1 must be quantized 

at 271 • 
at Nj . 


related Ginzburg-Landau action) with ^ f b 1 Ada 1 A 
ip f a 2 A a 1 A da 1 , where p = , k E Z n 12 . As a 

result, the total classification of irreducible SPT phases 
with Y[ 2 i -%-Nj symmetry is given by (Ztv 12 ) 2 - 

Physically, the quantization in Eq. (14) can be under¬ 
stood in the following way. By performing the gauge 
transformations (9), the total symmetry charges 

(/=l,2;/ = 3 — I) are not gauge invariant: 



Since the 1-form external gauge potential A 1 is quantized 
at we may enforce the additional terms are divisible 

by JVj, i.e. p f M3 d\ T A da 2 /Nj E Z such that the par¬ 
tition function is still invariant (e* 5 -A e lS+27Tl = e lS 
with the action S defined in Eq. (12)). Concretely 
speaking, p should be quantized properly in order to 
ensure that there exists an integer pair {rrA^m 2 ) such 
that the two equations below always hold: (/ = 1,2; 
1 = 3 — 1.) p f M 3 dx 1 A da 2 = Nimi . To proceed 
further, let us consider M 3 = S 1 x T 2 such that for 
1 = 1: 2tt p f M3 dx 2 A da 2 = 2tt p f sl dx 2 f T2 da 2 = 
2 ttp(27t) 2 x integers where Eqs. (8,10) are applied. There¬ 
fore, the existence of integer mi requires that “47r 2 p is 
divisable by Ni”. Likewise, we can obtain a similar con¬ 
dition: u 47r 2 p is divisable by N 2 ” . In summary, p is quan¬ 
tized as: p = ^2 N ^ 2 2 , k e Z , where ATi 2 is the greatest 
common divisor (GCD) of Ni and N 2 . 

However, k £ I doesn’t mean that the classification 
of bosonic SPT is necessarily Z. Let us consider the 
following shift operation : 


db 1 -A db 1 — — ^_^ 2 a 2 A da 2 , 
2tt7Vi2 

db 2 -A db 2 + — N}^ 2 a 1 A da 2 , 
2tt7Vi2 

k -A k + K 1 + K 2 , 


(16) 

(17) 

(18) 



= fc°’ ±v 


(13) 


D. Symmetry-protected quantization 

The implementation of the discrete symmetry leads to 
the following quantization and periodicity of p: 


where K 1 and K 2 are two integers parametrizing the 
above shift operation. Likewise, one can check that 
TQFT partition function is invariant under the above 
shift operation if and only if 


1 

27r 



after shift 


1 

27T 



before shift 


is divisible by Nj ( for / = 1, 2). 


(19) 


k NiN 2 
P ~A^ V 12 


k E Zat 12 , 


(14) 


where Ni 2 denotes the greatest common divisor (GCD) 
of Ni and N 2 . Eq. (14) indicates that there are totally 
7Vi2 topologically distinct SPT phases that are described 
by TQFT Eq. (6). For example, the SPT phase labeled 
by k is equivalent to that labeled by fc + Ni 2 . In addition 
to Eq. (6), we may also consider a TQFT (as well as its 


As a result, e lS -A e tS + 2 ' Kl = e lS with the action 
S defined in Eq. (12). By using the shift operations 
Eqs. (16,17), the two equations in Eq. (19) reduce to: 


K x NiN 2 

47T 2 7Vi2 

K 2 NiN 2 

47T 2 7Vi2 


/ 

/ 


a 2 A da 2 is divisible by N \, 
a 1 A da 2 is divisible by N 2 . 
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Integrating out b 1 leads to: f sl a 1 = — f sl A 1 = x 
0, dbl, • • • . By further using the Dirac conditions ( 8 ), we 
end up with two conditions: divisible by 

N\ and K A7T 2 ^ yyy is divisible by N 2 . The general solu¬ 
tion is K 1 = lNi 2 ,K 2 = VN 12 with /,/' E Z. Therefore, 
by using Bezout’s lemma, K 1 + K 2 is quantized at GCD 
of N 12 and JV 12 , he., N± 2 . In other words, the minimal 
shift of k [see Eq. (18)] is Ni 2 . For our bosonic SPT 
phases, let us start with k = 0 which describes the triv¬ 
ial phase. Then, by increasing k, distinct bosonic SPT 
phases labeled by k = 1 , 2 , • • • are obtained in succession. 
But once k = N\ 2 , the state returns to the trivial state 
immediately since it can be connected to k = 0 trivial 
state just by the above shift operation. One may also 
consider TQFT with action: 

2 

£tqft = 7 ^— ^ J b 1 Ada 1 A- ip ja 2 A 0 } A da 1 -\-Sm (20) 

which leads to another set of SPT states with p = 4 ^^ 
and classified by k E Z n 12 - As a result, the total number 
of distinct 3D SPT states with Y\ 2 Z^ T symmetry is given 
by the cyclic group (Zjv 12 ) 2 , which is also consistent to 
group cohomology: % 4 (n 2 I Nl , U(l)) = (l Nl2 ) 2 . 


IV. OTHER SYMMETRIES AND 
APPLICATIONS 


A. Other Abelian symmetries 


In Appendix A, the irreducible (see Table I) SPT 
phases with Z^ x x Zjv 2 x Zn 3 are derived, where the 
actions of TQFT are given by: 


. 3 


5tqft= /*A da Ij r ip Ja 1 A a 2 A da 3 + Sm , 

( 21 ) 


and, 


Stqft = — ^ Jb 1 A da Ij r ip Jc 


a 2 A a 6 A da 1 + Sm • 
( 22 ) 


( 21 ) and ( 22 ) produce (Zat 23 ) 2 classification of irreducible 
SPT phases. Since SPT phases with only two Zjv sym¬ 
metries can be viewed as reducible SPT phases with three 
Z;v symmetries, we conclude that from TQFT approach 
the full classification of SPT phases with G— IT, Z JVr 
symmetry is given by: 'H 4 (II? , U(l)) = (Zjv 12 ) 2 x 

(Zj v 23 ) 2 x (Zjv 13 ) 2 x (ZjVi 23 ) 2 • Therefore, we have achieved 
all SPT phases with either n? ° r n? Ztvj symmetries 
that were predicted in group cohomology. 


The irreducible SPT phases with r Ti Z Nl symmetry 
have been studied in Ref. 25 where TQFT is given by: 

5tqft = J b 1 Ada 1 + ip J a 1 Aa 2 Aa 3 A a 4 (23) 

with conditions Eq. ( 8 ) with / = 1, • • • ,4. p = 
with k E Zjv 1234 where V 1234 is GCD of all 
TV/’s. Thus, for arbitrary finite unitary Abelian group 
G = Z^Vi x Ztv 2 • • •, we can derive the TQFT descriptions 
for all SPT phases classified by the group cohomology: 


n\z Nl xz N2 --- ,u(i)) 

=n I<: j(Z NlJ ) 2 xU i<j<k (Z N ijk ) 2 xU i<j<k<l Z N ijkl . 

In fact, all the results for finite unitary Abelian group can 
be generalized into cases that involve U(l) symmetry as 
well. In Appendices, we show that there is no nontrivial 
SPT with U(l) /c or ZjVi x U(l) symmetry, and the ir¬ 
reducible SPT phase with x U(l) symmetry can 

be described by ^ f b 1 A da 1 + ip f a 1 A a 2 A da 3 , 
where p = ^ , k E Z Nl2 gives rise to a Z Nl2 

classification. Together with the (Ztv 12 ) 2 reducible SPT 
phases contributed from SPT phases protected merely 
by the ri/ZjVj part, we derive the TQFT descriptions 
for all the (Zn 12 ) 3 SPT phases predicted by group coho¬ 
mology classifications. The above scheme can be easily 
generalized to arbitrary Abelian group symmetry with 
the form G = Z Nl x Z N2 • • • x U(l) fc = G x U(l) fe , that 
is, in addition to the TQFT descriptions for finite uni¬ 
tary Abelian group piece G, one can always pick up a 
subgroup Zatj x Z x U(l) to construct Zn tj additional 
SPT phases. 


B. Connection to “decorated domain walls” 
picture 

Ref. 40 proposed a way called “decorated domain 
walls” to construct the ground state wave function of 
SPT states. Take Yl / Zjv 7 symmetry as an example, the 
SPT ground state can be formed as an equal weight su¬ 
perposition of 2D domain walls of Zn 1 on which a non¬ 
trivial 2D Ztv 2 SPT is placed. Here we use our TQFT 
approach to understand the “decorated domain walls” 
picture. Let us start with the action S in Eq. ( 12 ) and 
set A 1 = 0 , A 2 = A. We also add a dynamical Higgs 
term: 

Sm ggs = p-^BAdA, (24) 

where A and B are 1 -form and 2 -form dynamical gauge 
fields respectively. Sniggs breaks U(l) gauge symmetry 
of A 2 gauge field down to Zjv 2 - Integrating out B and 
b 2 leads to f M i a 2 = f M i A = x integer, a 1 A 

a 2 A da 2 becomes: ip f a 1 A A A dA. a 1 can be viewed 
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as domain wall of symmetry such that it may be 
replaced by where k\ G Z^q. In other words, given 
a Z^q symmetry domain wall “D(Z 7 \q)” labeled by Aq, 
there is a topological gauge theory on this 2D manifold: 


i r 27V2&1& 

Jv(z Nl ) N12 


A A dA , 


(25) 


where we have used p = 4-2 ^/v ^ 2 with /c G Zj\q 2 - By 
noting that 2N ^ k is even 18-24 , Eq. (25) suggests that 
the 2D state decorated on the domain wall is indeed a 
nontrivial 2 D SPT with Zj\r 2 symmetry. Such a way of 
thinking can be naturally extended to SPT phases with 
n? described by TQFT Eq. ( 6 ). For this purpose, we 
may introduce A and A that respectively gauge Z n 2 and 
Ztv 3 symmetry groups. As a result, on the Zj\q symmetry 
domain wall, we have the topological gauge action for 
Ztv 2 x Ztv 3 SPT phases: 


S[A,A] 


i 

47r 


L 


V(Z Nl ) 


2N 2 k 1 k 

W 2 


AA dA, 


(26) 


where Eq. (14) is applied and k G Zj\q 23 - 


C. Connection to 3 -loop statistics 

It was first shown in Ref. 62 that, after the global sym¬ 
metry Z ^q x Zw 2 x • • • of SPT phases is gauged at weak 
gauge coupling limit, the loop excitations of resultant 
gapped state exhibit the so-called 3-loop statistics. Based 
on our TQFT approach to SPT with G = Z^ xZjv 2 , we 
may formally gauge the symmetry by adding the coupling 
term: 

2 2 

■^coupling = J A dbI + flZ~^ BlAdAl ’ ( 27 ) 

where A 1 and B 1 are 1 -form and 2 -form dynamical gauge 
fields respectively. The first term is the minimal cou¬ 
pling term between A 1 and matter field current 
the second term is a Higgs term that breaks U(l) gauge 
symmetry of A 1 gauge field down to Z j\q • After integrat¬ 
ing out a 1 and b 1 , it turns out that the resultant new 
TQFT is given by (e.g. symmetry G= El? Ziv,): 

2 

S = Y2 J ^ B 1 A dA 1 +ipj A 1 A A 2 A dA 2 . (28) 

By integrating out B 1 , we will end up with an action for 
Dijkgraaf-Witten 74 gauge theory proposed in Ref. 50. We 
conjecture that Eq. (28) gives rise to the 3-loop statistics. 
Some progress has been made along this line 66-68 . 


V. CONCLUSIONS AND FUTURE 
DIRECTIONS 

In this paper, we attempt to provide a complete TQFT 
description for all 3D SPT phases with unitary Abelian 


group symmetry. Taking f | 7 Zj\q as an example, we start 
with a Ginzburg-Landau type action and illustrate the 
key mechanism via proliferating “nontrivial” symmetry 
domain walls for SPT phases with finite unitary Abelian 
group symmetry. We then rigorously derive the corre¬ 
sponding TQFT and compute the level quantization. Fi¬ 
nally, we consider generic 3D SPT phases with arbitrary 
unitary Abelian group symmetry, including U(l). All ir¬ 
reducible TQFT results for unitary Abelian group sym¬ 
metry are collected in Table I. In Appendices, we also dis¬ 
cuss possible TQFT descriptions with anti-unitary time 
reversal symmetry. Together with the previous work of 
bosonic topological insulators 55 , we provide a route to¬ 
wards a complete TQFT description for 3D SPT phases 
and reveal the key mechanism of these exotic quantum 
phases of matter. 

There are several interesting directions for future stud¬ 
ies. One of them is to construct the boundary of 3D SPT 
phases by using the present TQFT framework. There are 
two steps. First, we should consider TQFT in an open 
manifold and try to derive the boundary effective field 
theory. It is possible that the gauge invariant argument 
in boundary CFT derivation of Chern-Simons theory 2 is 
still applicable in the present TQFT. However, the ex¬ 
otic gauge transformations, e.g. Eq. (9), may lead to 
technical challenge. Second, the symmetry transforma¬ 
tions on the boundary effective field theory should reveal 
some nontrivial properties that forbid the realization of 
the boundary theory in a 2D plane alone (i.e. no 3D 
bulk). In addition, it is also interesting to study the 
symmetry-enriched topological phases (SET) by using 
TQFT. SET has topological excitations in the bulk and 
symmetry acts on excitations in a fractionalized man¬ 
ner (e.g. projective representation). For this purpose, 
one may replace ^ b 1 A da 1 by J ^ 7 A da 1 with 
ki = 2, 3, • • • and then study how to impose symmetry 
G. Along this line, one may study 3D SET phases with 
gauge group Z/ Cl x Z/~ 2 x • • • and symmetry group G. Fi¬ 
nally, it is interesting to extend all of the above studies to 
fermionic SPT and fermionic SET where the transparent 
quasiparticles are fermionic (e.g. electrons in FQHE). 
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Appendix A: TQFT of irreducible bosonic SPT 
phases with global symmetry Zn 1 X Zjv 2 x Z n 3 

In this Appendix, we derive TQFT and its classifica¬ 
tion when the symmetry is x Z n 2 x Z n 3 - The action 
is given by: 


Stqft = — f ^ ^ da T -\- ip j a 1 


A a 2 A da 3 + Sm 


(Al) 


The implementation of the discrete symmetry leads to 
the following quantization and periodicity of p : 


k N x N 2 
4^ n 12 


k G Ztv 123 i 


(A2) 


where N\ 2 (A123) denotes the greatest common divisor 
(GCD) of Ai and N 2 (A x , N 2 and A 3 ). Eq. (A2) in¬ 
dicates that there are totally A123 topologically distinct 
SPT phases that are described by TQFT Eq. (Al). For 
example, the SPT phase labeled by k is equivalent to that 
labeled by fc + A123. In addition to Eq. (Al), we may also 
consider a TQFT (as well as its related Ginzburg-Landau 
action) with ^ Y^i f b 1 A da 1 + ip f a 2 A a 3 A da 1 , where 
p = 4-2 N ^ 2 j k G Z n 123 • As a result, the total classifi¬ 
cation of irreducible SPT phases with Y\ 3 Zjv 7 symmetry 
is given by ( Z Nl23 ) 2 . 

Physically, the quantization in Eq. (A2) can be under¬ 
stood in the following way. By performing the gauge 
transformations ( 9 ), the total ZjVj symmetry charges 
(J=l, 2;/=3 — I) are not gauge invariant: 



Since the 1-form external gauge potential A 1 is quan¬ 
tized at we may enforce the additional terms are 

Aj x integer: p f M3 d\ T A da 3 = A/ x Z such that the 
partition function is still invariant. Concretely speak¬ 
ing, p should be quantized properly in order to ensure 
that there exists an integer pair (m^m 2 ) such that the 
two equations below always hold: (/ = 1 , 2 ; I = 3 — 1 .) 
p f M3 dx 1 A da 3 = Ni mi . To proceed further, let us 
consider M 3 = S 1 x T 2 such that for I = 1: 


27 vp 



A da 3 



= 27 rp( 27 r ) 2 x integers, 


where Eqs. ( 8 , 10 ) are applied. Therefore, the existence 
of integer mi requires that “47r 2 p is divisable by Ai” 
such that mi = x integer G Z. Likewise, we can 
obtain a similar condition: “47r 2 p is divisable by A 2 ”. 
In summary, p is quantized as: p = , fc G Z, 

where Ah2 is the greatest common divisor (GCD) of N\ 
and N 2 . 


However, fc G Z doesn’t mean that the classification 
of bosonic SPT is necessarily Z. Let us consider the 
following shift operation : 


,3 ,3 K 3 NiN2 i 2 

b 3 -> b 6 , T a 1 A a , 

2nN 12 

(A 4 ) 

db 1 — > db 1 — — A da 3 , 

2nN 12 

(A 5 ) 

db 2 ->• d6 2 x2 ^ liV2 da 3 A a i 
27 TAf 12 

(A6) 

k -> fc + a: 1 + a: 2 + a: 3 , 

(A 7 ) 


where K 1 , A 2 , and K 3 are three integers parametrizing 
the above shift operation. One can check that TQFT is 
formally invariant under the above shift operation. Fur¬ 
thermore, the shifts in b 1 fields should be consistent to 
the Z jsfj symmetry via the following relations: 


1 

27r 



after shift 



db 1 


before shift 


= Ni x integer ( for I = 1, 2, 3 ) (A8) 


which means that the change amount should be divis¬ 
ible by Ni for index I. By using the shift operations 
Eqs. (A 4 ,A 5 ,A 6 ), these relations reduce to: 


K 3 N iJV 2 
4tt 2 Ai2 
A 1 AiJV 2 
4tt 2 Ai2 
A 2 AiA 2 
4tt 2 Ai2 


J a 1 A da 2 + J da 1 A a 2 j is divisible by A 3 , 
J a 2 A da 3 is divisible by Ni , 

J da 3 A a 1 is divisible by N 2 . 


It means the integers (A 1 , A 2 , A 3 ) should be properly 
selected such that the three “integers” on the right hand 
sides always exist. Integrating out b 1 leads to: 



[ A 1 = A x 0, ±i, • • • . (A9) 

Js 1 A/ 


By using Eqs. ( 8 ,A 9 ), the above equations are reex¬ 
pressed as: 


K :i N-\ N-2 4-7T 2 (jV-| x integer + A r 2 x integer) 


47T 2 iVi2 

is divisible by N3 , 


NiN 2 


4 tt 2 . 


4.ir 2 N\2 N 2 


is divisible by Ni 


4 tt 2 . 


47T 2 A r i2 Ni 


is divisible by N2 ■ 


(A10) 

(All) 

(A12) 


By noting that (Ai x integer + A2 x integer) is always 
divisible by A i2 due to Bezout’s lemma, we obtain 
the following solution: A 3 /A 3 G Z, K 1 /N\ 2 G Z, 
and A 2 /Ai2 G Z. By using Bezout’s lemma again, 
we obtain the minimal shift of k is GCD(N\ 2 , A3) = 
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GCD(Ni, 7V 2 , AT 3 ) = 7Vi23. Thus, for our bosonic SPT 
phases, let us start with k = 0 which is trivial. Then, 
by increasing &, distinct bosonic SPT phases labeled 
by k = 1 , 2 ,-•• are obtained in succession. But once 
k = A/"i 23 , the state returns the trivial state immediately 
since it can be connected to k = 0 trivial state just by 
the above shift operation. 


Appendix B: TQFT of irreducible bosonic SPT 
phases with global symmetry {Zn 1 x • • •) x U(l) 

1. Zjv x U(l) 


Likewise, we can obtain the quantization of p: p = 
4^2 ? k eZ. We must also proceed further with the 

shift operation defined in Eqs. (A4,A5,A6,A7). By tak¬ 
ing account of the symmetry transformations Eq. (B4), 
Eq. (A 8 ) is changed to: 


w f dbl 

= Ni x integer( for 1 = 1 , 2 ) 


after shift 


1 

27T 



db 3 


after shift 


1 

27T 



before shift 

(B5) 


before shift 


= 0 
(B 6 ) 


In the following, we consider bosonic SPT phases with 
direct product of U(l) and several Zjv’s. First, let us 
consider G= Zn x U(l). In this case, we may consider 
TQFT with the conditions: 

Dirac conditions: — f da 1 £ Z, — f db 1 G Z 
2 ?r 2 tt 

(Bl) 

Symmetry: —- [ db 1 —>> —- [ db 1 + A (B 2 ) 
2 ^ Jm 3 2 tt 7^3 

which means that all bosons of / = 1 carry Z^v sym¬ 
metry while all bosons of / = 2 carry U(l) symme¬ 
try. The gauge transformations for both b 1 are given 
by: b 1 —>• b 1 — 2irp e IJ x J da 2 . Again, we require that the 
additional terms in b 1 fields after gauge transformations 
can be removed by the Z n symmetry transformations 
Eq. (B2). However, this can only be done for I = 1. The 
absence of I = 2 in Eq. (B2) leads to p = 0. 

In a similar way, we may show that bda + a 2 a 1 da 1 is 
also trivial. Note that bosonic SPT phases with either 
U(l) or Zjv symmetry are always trivial, we conclude that 
all bosonic SPT phases (both irreducible and reducible) 
with G= Zjv x U(l) symmetry are trivial. 


2 . n?^xU(l) 


Eq. (B 6 ) leads to: 

K 3 = 0 . (B7) 

By using the shift operations Eqs. (A5,A6), Eq. (B5) 
leads to 


1 2 J a 2 A da 3 is divisible by Ni , 
J da 3 A a 1 is divisible by N 2 . 


4tt 2 Ai2 

K 2 N x N 2 


4tt 2 Ai2 

By further using Eq. (A9) (/ = 1,2), we end up with: 


K 1 N\N 2 47t 2 . 


47T 2 A/i2 N 2 


is divisible by N\ 


K 1 N\N 2 47t 2 . . . U1 u Ar 

~ OAr -— is divisible by N 2 . 

4tt 2 Ai2 N x j 


which means that: K 1 /Ni 2 G Z,K 2 jN\ 2 G Z. Thus 
the minimal shift of k is Ni 2 . Therefore, the irreducible 
bosonic SPT phases are classified by cyclic group Zn 12 . 

In summary, the complete classification of bosonic SPT 
phases with f\ 2 j Z^ T x U(l) symmetry is given by Zn 12 x 
(Z tv 12 ) 2 = (Z 1 2 ) 3 where the additional two Z^ 12 parts 
arise from bosonic SPT phases with Y\ 2 Z symmetry 
only. 


Next, we consider G= Y[ 2 x U(l). Since there are 
three independent symmetry groups, we need to consider 
TQFT in Eq. ( 6 ). If we assume that the bosons of either 
/ = 1 or / = 2 carry U(l) symmetry, then, it is still 
concluded that no irreducible bosonic SPT phases exist. 
Therefore, it is sufficient to only explore the possibility 
of nontrivial irreducible bosonic SPT phases where the 
bosons of / = 3 carry U(l) symmetry. Then, the follow¬ 
ing conditions should be imposed: 


uirac 


Symmetry: 

( 1 = 1 , 2 ) 


Z7r JM 2 * 

c [ db 1 f i 
2tt Jm 3 2tt J 3 


Appendix C: TQFT of irreducible bosonic SPT 
phases with time-reversal symmetry (Z 2 ) 

1. Z 2 and U(l) x Z 2 

The complete classification of bosonic SPT phases with 
Z 2 (time-reversal symmetry with T 2 = 1 ) is given by 
(Z 2 ) 217 . One Z 2 index corresponds to bosonic SPT 
phases with surface “all-fermion topological order” 30 . 
Its TQFT description is proposed to be the form of 
u b A da + b A 6 ” in Ref. 55. The other Z 2 index corre¬ 
sponds to bosonic SPT phases with surface Z 2 topolog¬ 
ical order where both e and m carry Kramers’ doublet. 
Its TQFT description is just a multi-component b A da 
theory of level -1 but with unusual definition of bulk time- 
reversal symmetry transformation, as shown in Section 
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VII of Ref. 55. These bosonic SPT can be viewed as 
either bosonic topological insulators where U(l) symme¬ 
try doesn’t play role of symmetry protection or “bosonic 
topological superconductors” where U(l) symmetry is 
completely broken. 

Formally, U^xiZ^ is a non-Abelian symmetry group 
due to the semi-product operation “x” but we still sum¬ 
marize known results below. The complete classification 
of bosonic SPT phases with U(l)xZj symmetry is given 
by (Z 2 ) 317 , among which there are two Z 2 indices are 
given by bosonic SPT phases with merely time-reversal 
symmetry. The remaining one Z 2 index labels irreducible 
bosonic SPT phases with U(l)xZ^ symmetry so that 
both U(l) and Z^ play nontrivial role of symmetry pro¬ 
tection. In other words, both U(l) and Z^ symmetry 
play nontrivial role. It can be understood through Wit¬ 
ten effect discussed in Refs. 55 and 56. The TQFT de¬ 
scription of this Z 2 index is given by a multi-component 
b A da in Sec. VI of Ref. 55. 


2. U(l) x Tl 

The complete classification of bosonic SPT phases with 
U(l) x IT} symmetry is given by (Z 2 ) 417 , among which 
there are two Z 2 indices are given by bosonic SPT phases 
in C 1. The remaining two Z 2 indices label irreducible 
bosonic SPT phases with U(l) x Z^ symmetry. 

One Z 2 index can be understood in the similar way 
to the Witten effect in bosonic topological insulators 
discussed in Refs. 55 and 56. The response theory is 
still given by O = 2tt F A F response action where the 
external gauge field A here is “spin gauge field” that 
is pseudo-vector. Under time-reversal transformation 
“electric field” E changes sign while “magnetic field” B 
not 30 . The TQFT description of this Z 2 index is the 
same as the multi-component b A da in Ref. 55 by just 
changing the time-reversal transformation of external A 
from polar-like to pseudo-like transformations. 

The other Z 2 index is signalled by surface Z 2 topo¬ 
logical order where e carries half charge (compared to 
the fundamental U(l) charge unit of bulk bosons) and m 
carries Kramers’ doublet 54 . For the purpose of TQFT de¬ 
scription, we still start with the following two-component 
b A da theory with a coupling to external “spin gauge 
field” A: 

S = [ i—b 1 Ada 2 + i [ — b 2 A da 1 

J 2tt J 2tt 

+ i J — 6 2 A da 2 + ^coupling • (Cl) 

It corresponds to ^ l -K IJ Jb A da with K = (J^). 


^coupling is given by: (gi, g 2 e Z , F = dA ) 

^coupling = J l -^F A da 2 + ^b 2 A dA . (C2) 

Time-reversal transformation is defined as: (/ = 1,2; 
=x,y,z) 

Ta^l 1 = al ,Ta{ T” 1 = -a? , (C3) 

Tblj- 1 = -b^ THjT- 1 = b! 0 . (C4) 

Integrating out a 1 leads to local flatness of b 2 , which 
can be resolved by introducing a new 1-forma gauge field: 
b 2 = da 2 . Thus, the time-reversal transformation is given 
by: 

TdlT- 1 = -d 2 0 , Td{ T- 1 = d 2 . (C5) 

The term ~ b 2 A da 2 in Eq. (Cl) provides a surface 
Chern-Simons term: 

c d = i-e^ x d 2 d u a 2 (C6) 

7T 

which can also be reformulated by introducing a ma¬ 
trix Kq - e - (°o) i n the standard convention of K- 
matrix Chern-Simons theory. 2 a 2 and a 2 form a 2- 
dimensional vector (a 2 ,a 2 ) T . The ground state of 
Eq. (C6) supports a Z 2 topological order associated with 
four gapped quasiparticle excitations (1, e, m, e:). By us¬ 
ing b 2 = da 2 , Eq. (C2) reduces to its surface counterpart: 
| ^ x A,d v a\ + ge^A^al. 

Based on the Chern-Simons term in Eq. (C6), one 
may calculate the electric charge carried by each quasi- 
particle: Q e = (q 1 ,q 2 ){K 9 )~ 1 (l,0) T = f, Q m = 
(QijQ2){Kg)~ 1 (0,l) T = ■ Physically, both e and m 

quasiparticles can always attach trivial identity particles 
to change their charges by arbitrary integer so that q\ and 
#2 are integers mod 2, namely, q\ ~ q\ + 2 , q 2 ~ #2 + 2 . If 
we consider (# 1 ,^ 2 ) = (2,1), then Q e = \ , Q m = 1 indi¬ 
cating that e carries half-charge of U(l) symmetry. Since 
a 2 transforms as a pseudo-vector under time-reversal 
transformation defined in Eq. (C5), its gauge charge, m 
particle, is able to carry Kramers’ doublet. This sur¬ 
face is what we need: e carries half-charge while m car¬ 
ries Kramers’ doublet. We like to point out that the 
surface state cannot be realized on 2D plane alone un¬ 
less time-reversal symmetry is broken. More concretely, 
let us investigate the coupling term LL ^e^ vX A ll d v d\ that 
changes sign under time-reversal transformation since A 
is a pseudo-like vector: TAqT -1 = — Ao ,7" AiT -1 = Ai . 
Therefore, this state necessarily break time-reversal on 
2D plane. Since the sign change can be remedied by 
shifting q 2 through the identification “gi ^ qi + 2 , q 2 ~ 
q 2 + 2 ”, the surface state is time-reversal invariant. Note 
that is the external gauge field is the usual electromag¬ 
netic field that is a polar-like vector, then, the symmetry 
group is U(l) x Z 2 and there are no obstruction for re¬ 
alization of the above surface state on 2D plane. 
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